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A black universe is a nonsingular black hole where, beyond the horizon, there is an expanding, asymptotically isotropic 
universe. Such spherically symmetric configurations have been recently found as solutions to the Einstein equations 
with phantom scalar fields (with negative kinetic energy) as sources of gravity. They have a Schwarzschild-like causal 
structure but a de Sitter infinity instead of a singularity. It is attempted to obtain similar configurations without 
phantoms, in the framework of an RS2 type brane world scenario, considering the modified Einstein equations that 
describe gravity on the brane. By building an explicit example, it is shown that black-universe solutions can be 
obtained there in the presence of a scalar field with positive kinetic energy and a nonzero potential. 
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1. Introduction 



The problem of singularities is one of the long-standing 
problems in the classical theories of gravity. Singulari- 
ties are places where general relativity or another clas- 
sical theory of gravity does not work. Therefore, a full 
understanding of the physics of phenomena under study 
(origin and fate of our Universe, gravitational collapse 
etc.) requires avoidance of singularities or/and modifi- 
cation of the corresponding classical theory or address- 
ing quantum effects. There have been numerous at- 
tempts on this trend, some of them suggesting that sin- 
gularities inside the event horizons of black holes should 
be replaced with a kind of regular core ( [1], see [2] for 
a recent review), others describing bouncing or "emer- 
gent" universes (see, e.g., [3,4] for reviews). 

In our view, of particular interest are models which 
combine avoidance of singularities in both black holes 
and cosmology, those which have been termed black uni- 
verses [5, 6]. These are regular spherically symmet- 
ric black holes, with the same causal structure as the 
Schwarzschild black hole, but where a possible explorer, 
after crossing the event horizon, gets into an expand- 
ing universe instead of a singularity. Thus such hypo- 
thetic configurations combine the properties of a worm- 
hole (absence of a center, a regular minimum of the area 
function) and a black hole (a Killing horizon separating 
R and T regions). Moreover, the Kantowski-Sachs cos- 
mology in the T region is asymptotically isotropic and 
approaches a de Sitter regime of expansion, which makes 
such models potentially viable as models of our acceler- 
ating Universe. 
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Such objects are naturally obtained [5,6] if one con- 
siders local concentrations of dark energy (DE) repre- 
sented by different forms of phantom matter: phantom 
scalars in scalar-tensor theories of gravity and the so- 
called k-essence whose most general Lagrangians have 
the form F((f>,X) where X = V Q </>V a <^. In all such 
cases, there is at least one classical field with a negative 
kinetic term. 

The phantom behavior of DE (such that its pressure 
to density ratio w = p/e < —1) is favored by modern 
cosmological observations. More precisely, the "Gold" 
supernova sample data [7] slightly favor a phantom be- 
havior of DE at small redshifts z < 0.3 along with cross- 
ing the phantom divide w = — 1 at larger z [8] (see [9] 
for further references) . Other cosmological data suggest 
a cosmological constant as the best fit but still do not 
exclude recent phantom DE behavior, see [10]. The lat- 
est estimates of w also peak somewhat near —1 and 
admit w < — 1 [11]. 

Meanwhile, there exist models that admit a phan- 
tom DE behavior without explicitly introducing phan- 
tom fields. Among them the simplest is the generic 
scalar-tensor gravity with non-zero scalar field poten- 
tials (see, e.g., [12]) which has sufficient freedom to de- 
scribe all observational data. 

There are theoretical reasons for considering phan- 
tom fields: they naturally appear in some models of 
string theory [13], supergravities [14] and theories in 
more than 11 dimensions like F-theory [15]. 

Nevertheless, one should bear in mind that a classical 
field with a negative kinetic term can have an arbitrar- 
ily large negative energy of high-frequency oscillations, 
which is quite undesirable from the viewpoint of quan- 
tum field theory: it can lead to runaway production of 
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particles and antiparticles accompanied by production 
of equal negative energy of the phantom field itself (see, 
e.g., [16]). Nothing of this kind is observed, which casts 
serious doubt on possible existence of phantom fields. 
Moreover, as was recently argued in [12], cosmological 
models with a phantom scalar field cannot explain the 
observed large-scale homogeneity and isotropy of the 
Universe. 

Thus there exist arguments both pro et contra phan- 
tom fields, and the latter seem somewhat stronger. In 
any case, it is reasonable to try to avoid such fields in 
modelling real or hypothetic phenomena. 

Accordingly, we here try to show that black-universe 
models can be obtained without invoking phantom 
fields. This appears to be possible in the framework 
of the brane world scenario, using the modified Einstein 
equations [17] describing gravity on the branc. 

The brane world concept describes our world as a 
4D surface (brane) supporting all or almost all matter 
fields and embedded in a higher-dimensional space-time 
(called the bulk). This concept traces back to the early 
80s and leads to a variety of consequences in cosmology, 
gravitational and particle physics, see the reviews [18]. 
In particular, brane worlds turn out to be a natural 
framework for wormholes [19, 20] (see also references 
therein) since there the modified Einstein equations [17] 
[see Eqs. (13)] contain a source term of geometric 
origin which need not observe the usual energy condi- 
tions. And, as we shall see in the present paper, it is this 
source term that can replace phantom fields in building 
black-universe models. 

The modified Einstein equations (EE4) used here 
correspond to the so-called RS2 scenario: a single brane 
in a Z2 -symmetric 5-dimcnsional bulk, with all fields 
except gravity confined on the brane. It generalizes 
the second Randall-Sundrum model comprising a single 
Minkowski brane in an anti-de Sitter (AdS) bulk [21]. 
However, in other brane-world scenarios, the effective 
4D Einstein equations also contain terms similar to E v , 
e.g., on codimcnsion-1 branes without Zi symmetry [22] 
and in Dvali-Gabadadze-Porrati brane worlds [23] with 
different kinds of induced gravity terms [24]. Thus we 
can anticipate that black-universe solutions similar to 
ours exist in such brane worlds as well, though proba- 
bly under some other conditions. 

The paper is organized as follows. In Section 2 we 
briefly discuss the Einstein-scalar field equations for 
static, spherically symmetric systems in general rela- 
tivity and present the simplest black-universe solution 
with a phantom scalar. In Section 3 we analyze the 
similar set of equations in a brane world and show that 
for our purpose we can neglect some terms, namely, 
those quadratic with respect to the stress-energy tensor 
of matter. Section 4 is devoted to attempts to obtain 
black-universe solutions by properly employing the free- 



dom that exists in this system. Section 5 contains some 
concluding remarks. 

2. Black-universe solutions in general 
relativity 

Consider the general static, spherically symmetric met- 
ric 



ds 2 =A(u)dt 2 - ^ " 



A{u) 



r 2 {u)dVt 2 



(1) 



where d£l 2 = d9 2 + sin 2 ddtp 2 is the metric on a unit 
sphere. 3 The metric (1) is written in terms of the 
"quasiglobal" coordinate u , which is particularly conve- 
nient for dealing with Killing horizons where it behaves 
in the same way as the manifestly well-behaved Kruskal- 
likc null coordinates. For this reason, in terms of u, one 
may consider regions on both sides of such a horizon 
remaining in a formally static framework. 

The two functions, A(u) (often called the redshift 
function) and r{u) (the area function, equal to the ra- 
dius of a coordinate sphere at given u) entirely deter- 
mine the geometry under consideration. Horizons cor- 
respond to regular zeros of A(u) . 

Our interest is to find black-universe solutions which, 
by definition, must have the following properties: 

1. Regularity in the whole range u € R; 

2. Asymptotic flatness as u — > — 00 (without loss of 
generality), i.e., r(u) w —u, A(u) 

3. A de Sitter asymptotic as in 
region (A < 0) where r(u) <~ u, — A(u) 

4. A single simple horizon (i.e., a simple zero of 
A(u) ) at finite u. 

As shown in [6], such solutions can exist in general 
relativity with various phantom sources. Let us here 
present the simplest example with a minimally coupled 
scalar field having the Lagrangian 



-> i; 

+00, i.e., a T 
2 . 



L^-ed a ^d a <j>~V(<j>), 



(2) 



where e = +1 corresponds to normal scalar fields with 
positive kinetic energy, s = —1 to phantom fields, and 
V{4>) is a potential. The set of Einstein-scalar equations 
for the metric (1) and <f> = 4>{u) may be written in the 
form 



e(ArV)' = 

(A'r 2 )' = 

2r"/r = 

A(r 2 )"-r 2 A" = 



-r z dV/dcj), 
-2r 2 V; 

J/2. 



(3) 
(4) 
(5) 
(6) 



3 Our sign conventions are as follows: the metric signature (H — 

) ; the curvature tensor R" ' Mp „ = f „r^ p — . . . , so that, e.g., the 

Ricci scalar R > for de Sitter space-time, and the stress-energy 
tensor (SET) such that T t * is the energy density. 
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where the prime denotes dj du . The scalar field equation 
(3) is a consequence of Eqs. (23)-(25), which, given a 
potential V(4>), form a determined set of equations for 
the unknowns r(u), A(u), <j>(u). 

As is evident from (24), black- universe solutions can- 
not be obtained with e — +1 because in this case r" < 
which is incompatible with requirements 1-3 (instead of 
u G M, there will be inevitably r = at some finite 
u, which is either a singularity or, at best, a regular 
center). 

A particular solution to these equations with e = — 1 
is given by [5] 



B(u) 



r = (u 2 + 6 2 ) 1/2 , b = const > 0; 

A(u) _ c 1 
r 2 (u) ~ b 2 b 2 + u 2 

u ( bu u\ 

±\/2 arctan(w/6) + 0o, 

c r 2 + 2u 2 
b 2 r 2 

u ( 3bu r 2 + 2u 2 

£j3 I ljr»2, rp2. 



v = - 



U 

arctan — 



(7) 

(8) 
(9) 

(10) 



with c, uq, 4>o = const. In particular, 



B{±oo) = --V{±oo) = —^—. 



(11) 



Choosing in (9), without loss of generality, the plus sign 
and 0o = 0, we obtain for V{4>) 



Uq 



— ^3 [3 sin ip cos ip + -0(3 — 2 cos 2 0)] , 



V2' 



(12) 



The solution is asymptotically flat at large negative 
u ( B — > , A — > 1 ) under the condition 26c = ttuq , 
and the Schwarzschild mass, defined in the usual way, is 
then m = — uo/3. Then, assuming m > 0, we find that 
B(+oo) = — 'Sirm/b 3 = const < 0, which corresponds 
to a de Sitter asymptotic with a cosmological constant 
A > 0. The horizon position is found by solving the 
transcendental equation B(u) = 0. 

This is an example of a black-universe solution. 
Other cases of the solution (7)-(10) include wormholcs 
and asymptotically AdS configurations, see more details 
in [5,6]. 



3. Gravity on the brane 

3.1. Modified Einstein equations 

The gravitational field on the brane is described by the 
modified Einstein equations [17] 



4K 



(13) 



where G V U =R V ,- h5"R is the 4D Einstein tensor, A 4 is 
the 4D cosmological constant expressed in terms of A5 
and the brane tension A : 



1 



1 



(14) 



A 4 = - (^A 5 + -xgA 

x| = 8ttGn = x|A/(67r) = mj is the 4D gravitational 
constant; Gat is Newton's constant of gravity, and 7714 
is the 4D Planck mass; 

k 5 = m 5 3 ^ 2 , m 5 being the 5D Planck energy scale; 

is the SET of matter trapped on the brane; 

11^ is a tensor quadratic in , obtained from matching 
the 5D metric across the brane: 

K = i T » T « \ TT l - ¥1 (, T ^ Ta " \ T2 ) - (is) 

where T = ; lastly, E» is the "electric" part of the 
5D Weyl tensor projected onto the brane: in proper 5D 
coordinates, 

E»u = 5^^C ABC Dn B n D (16) 

where A,B,... are 5D indices and n A is the unit normal 
to the brane. By construction, is traceless, £^ 



0. 

Other characteristics of Ef. are unknown without 
specifying the 5D metric, hence the set of equations (13) 
is not closed. In isotropic cosmology this leads to an ad- 
ditional arbitrary constant in the field equations, con- 
nected with the density of "dark radiation" [18]. For 
static, spherically symmetric systems to be discussed 
here, this freedom is expressed in the existence of one 
arbitrary function of the radial coordinate. 

3.2. Reasons for neglecting 11^ 

Let us show that under quite reasonable conditions we 
can neglect the tensor IPJ in (13). 
We put A 4 = 0, so that 



I A 5 



>%A = 6-7T (X4/X5) 



(17) 



and use the observational restriction on the bulk length 
scale t which follows from the recent short-range New- 
tonian gravity tests [25], showing that Newton's inverse- 
square law hold at length scales greater than about 0.1 
mm. This means that if we live on an RS2-likc brane, 
the bulk length scale can be estimated as 



£=(6/|A 5 |) 1 /2<io- 2 cm. 



(18) 
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Note that the 4D Planck scale in our notation is 

™ 4 = *l 1 = (8ttG n « 2.4-10 18 GcV, 
Z4 = l/m,4 = i=a 8T0 cm. 

Combining (17) and (18), we obtain 

m 5 /m 4 = (7r£/Z 4 ) _1/3 > 10 -10 , (19) 

so that the 5D Planck energy scale in this scenario is at 
least about 10 8 GeV. 

Now, we can assert that the term with 11^ is neg- 
ligible in (13) as compared with the term as long 
as 

k\W 2 < k\W, 

(where W characterizes the magnitude of T£ , say, the 
absolute value of the largest component of T£ ) , or 

W < mjj/m| = mt(m 5 /m 4 ) 6 , (20) 

where m\ w 3.5-10 73 GeV 4 w 8.4-10 90 g • cm" 3 is the 
Planck density while the second factor is, according to 
the experimental bound (19), about 10~ 60 or larger. As 
a result, for the "density" W we have 

W < 10 30 g • cm" 3 . 

Recalling that the density of nuclear matter is about 
10 13 g • cm -3 , it is clear that this bound certainly holds 
for any thinkable matter. 

3.3. Brane gravity with a scalar field 

Consider Eqs. (13) for static, spherically symmetric con- 
figurations of a normal scalar field with an arbitrary 
potential, neglecting the term IF . So the scalar field 
Lagrangian has the form (2) with e = 1 . The general 
static, spherically symmetric metric is again taken in 
the form (1), with the quasiglobal coordinate u. 

The scalar field EMT is conservative, so the same is 
required for E£ . If we take it, for convenience, in the 
form 

E; = diag(-P-A/, -P, P+Af/2, P+Af/2), 

(21) 

where P and / are some functions of the radial coordi- 
nate u (so that, as required, its trace is zero), then the 
conservation law V a E" is written as 

(Pr 4 )' = \r 6 f{A/r 2 )'. (22) 

Eqs. (13) may be written in the form 

^(A'r 2 )' = -V - P - Af; (23) 
2r"/r = - f + /; (24) 
A(r 2 )" - r 2 A" - 2 = 2P + ^Af. (25) 



The scalar field equation (Ar 2 4>')' — r 2 dV/d(j> follows 
from (23)-(25) combined with (22). 

Thus, if V{4>) is specified, we have four indepen- 
dent equations (22)-(25) for five unknown functions of 
u: <p, A, r, f and P. The system becomes still more 
undcrdctcrmincd if V{<p) is not specified: in this case we 
can choose as many as two functions by hand to obtain 
a solution. 

4. Attempts to obtain black-universe 
solutions 

4.1. Models with r" > 

To have a regular positive minimum of r{u) , we must 
have r" > at least in some range of u. By (24), this 
can be achieved only with / > . We first try to obtain 
such a solution with r" > everywhere, i.e., / > cf)' 2 . 
To integrate (22), let us suppose 

/ = 2C/r 6 , C = const. (26) 

Then we have 

Pr 4 = CB + C\, B(u) = A/r 2 , 

d = const. (27) 

Putting for simplicity C\ = and substituting (27) into 
(25), we obtain an equation connecting B(u) and r{u) : 

{r A B')' + 2 + IQCB/r 2 = 0. (28) 

If both r and B are known, the quantities (f> and V are 
found from (23) and (24). Thus it remains to choose r 
and B satisfying (28) and providing the above proper- 
ties 1-4. 

This task turns out to be hard if at all solvable. 
Thus, if we choose r 2 = u 2 + b 2 (b = const > 0) with 
good asymptotics at large u, we obtain r"/r — b 2 /r 4 , 
whence (24) combined with (26) leads to <p' 2 < at 
large u, i.e., a contradiction. 

Another choice, r 4 = u 4 +6 4 , is compatible with (24) 
if C > 3b 4 . However, a numerical solution of Eq. (28) 
leads to functions B{u) with at least 3 zeros, in other 
words, at least 3 horizons in space-time (see an example 
in Fig. 1). 

The nature of Eq. (28) apparently implies an oscilla- 
tory behavior of B(u) if the coefficient 10C/r 2 (origi- 
nating from E^ ) is large enough in a sufficiently wide 
range of u . A possible way of avoiding such a behavior 
is to try a "more concentrated" distribution of / and 
P. 

4.2. Models with / - S(u) 

Let us try to find a black-universe solution assuming 
that f(u) is concentrated on a single sphere, e.g., u = 0. 
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Figure 1: Example of an oscillating solution for B(x) 

We now make the equations dimensionless by putting 
x = u/b, Pr 4 = Q{x)b 2 , fr 4 = 2F(x)b 2 1 



r = br(x), B = B(x)/b 2 



(29) 



where b = const > specifies a length scale of the con- 
figuration, and to justify our rejection of 11^ we assume 
b ^S> I (see above) . 

In what follows we omit the bars over r and B. 
Eqs. (22), (24) and (25) take the form 



Q' = r 2 FB', 

r" l/2 F 
— = + - 
r 

(r 4 B')' + 2+^- 



6BF = 0, 



(30) 
(31) 
(32) 



where the prime denotes d/dx and ip = <j>/y/2. 

To begin with, we choose r(x) so that r" < at all 
x ^ 0: 



r 2 {x) = (|z|+c) 2 -l, 



c = const > 1. 



(33) 



This conforms to both flat and de Sitter asymptotics. 
Then, in Eq. (31), the quantity r"/r has a delta-like 
singularity at x = . To compensate it and make rp' 
continuous at x = , we put 



F(x) = 2cr 2 5(x), 
We have then 
^i 1 



'"o 



r 2 (\x\ + c) 2 — 1 ' 



(34) 



(35) 



whence, without loss of generality. 
h(x) := — In 



h(c — x), x < 0, 

2h(c) -h(c + x), x>0, 



1 



x + 1 



X- 1 



(36) 



B(x) 
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Figure 2: The function B(x) obtained from a delta- like 
/(«) 



Eq. (30) gives for Q(x) : 

Q{x) = Q + 2crtB x0 6{x), 
where Qq = const, B x q — B'(0) and 



(37) 



6{x) 



1, x > 0, 
0, x < 0. 



is the Heaviside function. 

It remains to solve Eq. (32) and to find V . Eq. (32) 
now has the form 



(38) 



where Bo — B(0). We solve (38), choosing the initial 
conditions at large negative x in the Schwarzschild form 

2m 

A(x) = r 2 (x)B(x) = 1 - — , m = const. (39) 

\x\ 

Integrating (38) once, we obtain 

r 4 B' = 2{x - x) - 4Q h(c - x) 
- 89(x)cr 4 B xa [h{c) - h{x+c)} - l2B a cr 2 6(x), (40) 

where x$ = — m/3, see (39). The constant B x q is found 

as 



TqB x0 = 2x Q - 4:Q h(c) - 6B cr 2 , 



(41) 



where we have used the principal value of 9(0) equal 
to 1/2. Note that the term with B x q in (40) does not 



(i 




Figure 3: The area function r(x) chosen with a delta- 
like (a) and smooth (b) function f(r) 



contribute to the expression (41) because it vanishes at 
x = 0. 

To find B{x) , we first integrate (40) from — oo to 
x < 0, obtaining 



B(x) 



= -5 + (xo-c) 



c — X 



r 2 (x) 



h(c — x) 



1 - 2(c - x)h(c - x) 
r 2 {x) 



+ h 2 (c-x) 



(42) 



This also gives the constant Bq = B(0) , 
1 



Bn 



1 - 2cft(c) 



h(c) 



h 2 (c) 



(43) 



which may be used to obtain B{x) at positive x from 
(40) by integration from to x > : 



B(x) 



1 1 

Bq H 9 ^ 



x>o Tq 
+ [(so + c) - 4/i(c)(Qo + cr^o) - 6S oC r2] 
c + a; 



+ (Q + 2cr^ a o) 



1 - 2(c + i)/i(c + a;) 



r 2 (cc) 



+ h 2 (c+x) 



(44) 



where [/(a;)]* := f(b) — /(a). This solution (see Fig. 
2) really describes a black universe, but the delta-like 
distribution of the effective exotic matter, characterized 
by f(u) , causes an undesirable discontinuity of B' at 
1 = 0. 

The expression for the potential V(x) is rather cum- 
bersome and will not be presented here. 

4.3. A model with smoothed f(u) 

Evidently the qualitative behavior of the model will not 
change if we replace the delta-like distribution of f(x) 



B(x) 
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c=0.8;d=2;Q0=-1.2;x0=-2 



Figure 4: The function B(x) obtained from a smoothed 
but peaked f(u) 



with a smooth one but sufficiently peaked near x = 0. 
We will present an example of such a solution. Namely, 
let us preserve the notations (29), so that the field equa- 
tions have the form (30)-(32); however, instead of (33), 
we choose the following function r(x) : 



r 2 {x) = (\x\ + l) 2 -1, \x\>d, 



r 2 (x) = ax 2 + s 



d+1 



\x\ < d, 



(45) 
(46) 



where d — const > is sufficiently small and the con- 
stants in (46) are chosen to make r and r' continuous 
at x = ±d, as shown in Fig. 3. 

Let us take f(x) = at |x| > d and, for \x\ < d, as 
in (26), / = 2C/r 6 . Then from (30) we find Q(x) as 
follows: 




[B(x) - B(~d)}C, 
[B(d) - B(-d))C, 



x < —d; 
\x\ < d; 
x > d, 



(47) 



where Qo is an integration constant. The constant 
C (more precisely, its dimensionless counterpart C — 
C/b 4 -) is determined from Eq. (31), where we require 
continuity of ip' at x = ±d. Eq. (02') for the function 
B(x) = A/r 2 is then solved analytically for |x| > d but 
only numerically for |x| < d. 

As a result, at x < —d we obtain B(x) in the form 
(42) with c — 1. At x > d, an expression for B{x) is 
similar to (44) but obtained with initial conditions for 
B and B' at x = d that follow from numerical inte- 
gration of Eq. (32) in the range |x| < d. The results 
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Figure 5: The potential V(x) obtained from a smoothed 
but peaked f(u) 



are shown in Fig. 4. Fig. 5 shows the corresponding po- 
tential V(x) . Clearly these are black-universe models, 
where the asymptotic behavior of V(x) approaching a 
positive constant as x — + oo corresponds de Sitter ex- 
pansion with a positive cosmological constant. 

5. Concluding remarks 

We have built a family of black-universe solutions to the 
modified Einstein equations valid in an RS2 type brane 
world. They have been obtained without explicitly in- 
troducing any phantom matter. Just as was the case 
with wormhole solutions [19,26], the role of exotic mat- 
ter in the field equations is played by the "tidal" term 
of geometric origin, which has no reason to respect the 
energy conditions known for physically plausible matter 
fields. 

This new kind of solutions, having a black hole na- 
ture as seen from large negative x , supplements the sets 
of known examples of both spherically symmetric branc- 
world black holes (see, e.g., [27,28]) and black holes with 
scalar"hair" (see, e.g., [5,6,29]). 

Let us recall that the existence of black-universe 
models suggests the idea that our Universe could ap- 
pear from another, "mother" universe and undergo 
isotropization (e.g., due to particle creation) soon af- 
ter crossing the horizon. The Kantowski-Sachs nature 
of our Universe, as opposed to the more popular spa- 
tially flat models, is not excluded observationally [31] 
if its isotropization had happened early enough, before 



the last scattering epoch (at redshifts z > 1000). A 
tentative estimate obtained for cosmological evolution 
beginning from a horizon (the so-called Null Big Bang) 
in another kind of model, that with a vacuumlike static 
core, has shown [30] that such models isotropize very 
quickly and can now be quite observationally indistin- 
guishable from isotropic ones. One can also notice that 
we are thus facing one more mechanism of universes 
multiplication, in addition to the well-known mecha- 
nism existing in the chaotic inflation scenario. 

The presently obtained models do not pretend to 
be quite realistic, they simply show the possibility of 
such a scenario in principle. As any solutions to the 
effective 4D equations describing gravity on the brane, 
they certainly need extension to the bulk, whose finding 
is quite a challenging problem, although an extension 
always exists due to the known embedding theorems. 
Another problem to be solved is that of stability of such 
solutions. 
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